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A simple model of the two-state ratchet type is proposed for molecular chemical engines 
that convert chemical free energy into mechanical work and vice versa. The engine works by 
catalyzing a chemical reaction and turning a rotor. Analytical expressions are obtained for 
the dependences of rotation and reaction rates on the concentrations of reactant and product 
molecules, from which the performance of the engine is analyzed. In particular, the efficiency 
of energy transduction is discussed in some detail. 
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The Fi motor (Fi-ATPase), which is one of biological 
molecular motors, 1,2 is a remarkable molecular machine. 
It works as a rotary motor when it catalyzes the hy- 
drolysis of adenosine triphosphate (ATP) into adenosine 
diphosphate (ADP) and inorganic phosphate (Pi). 3,4 It 
can also generate (synthesize) ATP from ADP and Pi 
when its rotor is forced to rotate in the opposite direc- 
tion; 5,6 this is analogous to a heat engine working as 
a heat pump. A molecular machine, like the Fi motor, 
that can convert chemical energy into mechanical work 
and vice versa will be referred to as a molecular chemi- 
cal engine; it operates as a motor if it produces motion 
out of chemical energy, whereas it works as a generator 
if it generates "fuel" molecules of high chemical poten- 
tial from "waste" molecules of low chemical potential by 
consuming mechanical energy. 

It has been recognized that certain fundamental fea- 
tures of biological molecular motors can be captured by 
"Brownian motor" or ratchet models in which the sys- 
tem undergoes Brownian motion on a potential surface 
that changes stochastically between two or more profiles 
corresponding to different chemical states; see, for ex- 
ample, refs. 7-11 for reviews. For example, simple two- 
state ratchet models have demonstrated how the Brow- 
nian motion can be rectified to produce directed mo- 
tion; 12, 13 the dependence of the motor velocity on the 
concentration of the fuel molecule ATP have been an- 
alyzed with simple models 14-16 and with an elaborate 
model; 17 and the efficiency of energy transduction 18 and 
other measures of efficiency 19-21 have been discussed. Al- 
though the effects of the concentration of fuel molecule 
or the input free energy on the performance of molecu- 
lar motors have been studied in previous investigations, 
little attention has been payed to the effects of waste 
molecules. 

In this letter we propose a simple model for molecular 
chemical engines that explicitly takes into account the 
effects of both the reactant and product molecules. Our 
model is a variant of two-state ratchet models 7,8, 11-13 in 
general, and is closely related with the one outlined by 
Astumian and Bier 22 in particular. The main difference 
between our model and other two-state models worked 
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out earlier lies in the physics associated with transitions 
between the states: in the present model a transition oc- 
curs when the engine binds or releases a fuel or waste 
molecule, whereas it occurs as a result of ATP hydrolysis 
or the reverse reaction in some versions 14-16, 18 of previ- 
ous models; in another 22 the transitions represents bind- 
ing and unbinding of nucleotides, but the way how the 
nucleotide concentrations are related with the transition 
rates arc different from ours. Analytic expressions will 
be obtained for the dependences of rotation and reaction 
rates of the engine on the concentrations of reactant and 
product molecules. By using these results, the efficiency 
of energy transduction is studied for both the cases of 
the engine operating as a motor and as a generator. 

Our model engine has a rotor, whose rotational angle is 
denoted by 9. The conformation (structure) of the engine 
is assumed to change as 9 is varied. The engine catalyzes 
a hypothetical chemical reaction B instead of more 
complicated reaction ATP ADP + Pi taking place in 
biological molecular motors. The chemical potentials \xa 
and fj,B of molecules A and B, respectively, in the isother- 
mal environment of temperature T are assumed to satisfy 
the condition A/j, = fiA—fJ-B > 0; i.e., A is the fuel and B 
is the waste. The engine can bind at most one molecule 
A or B. For simplicity, the binding and dissociation of A 
or B are supposed to be possible only if is at a particu- 
lar value 9a or 9b, respectively, with < 9a < 9b < 2tt, 
i.e., if the engine is in a particular conformation. 

We adopt the following assumption 23 on the reaction 
taking place in the engine: the relaxation of the "reac- 
tion coordinate" describing this reaction is so quick that 
this coordinate is always in thermal equilibrium, and the 
forward (backward) reaction A — ► B (B — > A) proceeds 
with certainty when 9 is varied from 9a to 9b {9b to 9 a)- 
This will allow us to define an effective potential V\{9) 
(it is actually a free energy 23,24 ) of the engine when it 
carries a "ligand molecule," which is A for 9 ~ 9a and B 
for 9 ~ 9b- (Alternatively, the effective potential V\ may 
be reduced from more detailed models 17, 25,26 of the Fi 
motor, for example.) 

The engine is said to be in "state 1" or "state 0" de- 
pending, respectively, on whether it is occupied by a lig- 
and molecule or not. Let Vq(9) be the potential (free 
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Fig. 1. Schematic representation of the present model for a 
molecular chemical engine: an example of a pair of potentials 
Vq{6) and V\{ff) in states and 1, respectively, is shown together 
with possible pathways, indicated by arrows, which the engine 
can take. See the text for details. 

energy) of the engine in state 0. Then the torque exerted 
on the rotor by the engine in state j (j = 0, 1) is given by 
— dVj/dO. In Fig. 1, an example of the pair of potentials 
Vo and V\ (which is partly motivated by the analysis 27 
of experiments on the Fi motor) is shown, together with 
"pathways" (indicated by arrows) the engine can take. 
For simplicity, we assume that a large barrier of poten- 
tial V\ lies in the interval 9b<9<9 a + 2ti (mod 2tt) 
which cannot be surmounted by the rotor; we also as- 
sume that potential Vo has a barrier, which may be large 
or small, in the interval 9 a < 9 < 6b- If the latter bar- 
rier is so large that the dashed passes in Fig. 1 cannot 
be taken, a single forward revolution (increase in 9 by 
27r) of the rotor is always accompanied by a single for- 
ward chemical reaction A — > B; if this is the case, it is 
said that the "mechanochemical coupling" of the engine 
is tight. Otherwise, it is loose. 9 ' 28 

A transition between states and 1 occurs when the 
engine binds or releases a ligand molecule (vertical ar- 
rows in Fig. 1). Let Wj(6) be the rate of transition from 
state j to the other state at angle 9 of the rotor. The 
assumption mentioned above that the binding and the 
dissociation occur at particular values of 9 may be ex- 
pressed as 

Wj (e)=u;f5(6-e A )+wf6(6-e B ) (j = 0,l), (1) 

where and ui^ are positive constants associated with 
binding (j = 0) and dissociation (j = 1) of molecules 
A and B (see Fig. 1), and 5(9) is the delta func- 
tion. The use of the delta function is a mathematical 
idealization, which has been introduced by several au- 
thors 8, 14,15,19,29-31 to carry out various calculations an- 
alytically. 

The condition of detailed balance requires 24 the first 
equality in 

a = exp T^r — )= — = P<* (a = A,B) (2) 

\ k B T ) n u a 

to hold, where a = A or B, AV a = Vx(9 a ) - V (9 a ), 
and /cb is the Boltzmann constant. The second equality 
in eq. (2) comes from the fact that the concentration n a 
of molecule a is proportional to exp(p a /kBT) (n° is the 



value of n a for which \i a = AV a ), and the third equality 
defines the reduced concentration p a . The binding rate 
of molecule a will be assumed to be proportional to p a ; 
this implies that the dissociation rate is independent of 
p a due to eq. (2). Hence we have 

w « = n a p a , ujf = n a (a = A, B) (3) 

with some positive constant K a . 

We shall consider a situation in which a constant 
torque L > 0, opposing the forward rotation, is exerted 
on the rotor externally. If the rotor revolves forward, the 
engine does mechanical work against the load torque (it 
works as a motor). On the other hand, if it revolves back- 
ward and the backward reaction B — > A takes place, the 
engine works as a generator. 

The rotational motion of the rotor is supposed to 
be described by the Langevin equation in the over- 
damped limit. 10 Then, the rotation rate v [the average 
of (d0/d£)/27r] and the reaction rate r (the difference of 
the average numbers of forward and backward reactions 
per unit time) may be calculated from the steady-state 
solution to the Smoluchowski (Fokker-Planck) equation 

BP 87 

~jf • "-jI'j • "•: ,/': j (i = 0,l), (4) 

associated with the Langevin equation, for the prob- 
ability densities Pj(9,t) of 9 in state j at time t. In 
eq. (4), Wj is the transition rate given in eq. (1) and 
J-j = -D (dP j /d9) - {dVj/d9 + L)(Pj/^) is the prob- 
ability current density in state j with Do = k-^T/^ the 
diffusion coefficient associated with the free rotation of 
the rotor, and 7 the drag coefficient of the rotor. Pro- 
vided that the steady-state solution Pj{9) = Pj(9 + 2n) 
is normalized in such a way that J Q {P + Pi)d9 = 1, 
we have v = Jo + J\ and r = uoqPo{9a) — ^iPi(9a) = 
u^Pi(9b) — cJq Pq(0b)', note that Jo + J\ is independent 
of 9 in the steady state. 

The steady-state solution can be obtained analyti- 
cally 8, 14, 15, 19,29-31 for the Smoluchowski equation (4) 
with transition rates wj(9) given as sums of the delta 
functions. In the present model, in which the rate con- 
stants are given by eq. (3), the following results for v and 
r are obtained: 

v pa- (Jp B - a 
— = , (5) 

D Co + CAP A + CBPB + CABPAPB 

r _ b A p A - ab B PB 

Do C + C A p A + C B PB + CABPAPB ' 

where a = exp[(27ri — AV A + Ay^/feeT], and a, 6's and 
c's are dimensionless functions of load torque L (a, Co, ca 
and cb depend also on ha/ Do and kb/Dq). Coefficients 
a, &'s and c's are given by expressions, too complicated 
to be presented here, containing integrals of certain func- 
tions of Vj(9) + L9. A few properties of these coefficients 
should be mentioned: (i) &'s and c's are positive; (ii) a = 
and 6a = 6b for L = 0; (iii) a > and b A > bs for L > 0: 
and (iv) a — ► 0, b A — > 1, and bs — > 1 in the limit of large 
barrier of potential Vo ■ From the first two properties to- 
gether with the definitions (2) of p A and pb, we observe 
that v > and r > for Ap > in the absence of load 
torque (L = 0), as one would expect. 
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Fig. 2. The dependence of rotation rate v (solid line) and reaction 
rate r (dashed line) on load torque L for Ap/k^T = 15, 20, and 
25 with pb = 1.0. The parameters characterizing the potentials 
(Vb and Vi) and the transition rates are chosen as follows: Wo = 
Wi = Wk B T, <5> = 18k B T, $1 = 40fc B T\ 6 A = tt/3, 9 b = 
47r/3, 9c = 7T, and ka = k-b = Do- The upper inset is the 
magnification of a region where v ~ and r ~ for Ap/k B T = 
20. The lower inset shows the efficiency r) for Ap/k B T = 20. 



It is emphasized that the dependences of v and r on pa 
and given by eqs. (5) and (6) are quite general in that 
they are independent of potential profiles (Vb and VI), 
which affect only the values of coefficients (a, a, 6's and 
c's). We also note that both v and r depend on pa and ps 
separately while in the two-state models for molecular 
motors worked out previously 14-16, 18, 22 (and a single- 
state model 32 proposed recently) the chemical potentials 
come into play only through the difference Ap = pa~Pb 
(or Ap = patp — Madp — MPi if the reaction ATP <^ 
ADP + Pi is considered instead of A B). 

In what follows we discuss the properties of the engine 
extracted from eqs. (5) and (6) for a particular set of 
potentials Vb and V\. Here, for simplicity, we consider the 
piecewise linear functions shown in Fig. 1: the vertices of 
Vb are located at 6 = 9a, 9c, and 9 b with 9c satisfying 
9a < 9c < 9b, and those of V\ at 6 = 0, 9a, and 9b', 
the potential shapes are specified by parameters <f>o = 
V (9 C ) - V (9 A ), W = V (9 B ) - V (9 A ), *i = Vi(0) - 
ViiOs), and W 1 = VI (0 A ) - V^ B ). 33 

It can be shown that eqs. (5) and (6) are approximated 

by 



v w r « D pa/(c + c A pA + cabPaPb) (7) 

in the absence of the load (L = 0) if exp (Wj/ksT) 3> 1 
(j = 0, 1), exp [(Wo - <f> )/k B T] > 1, ka/D is not too 
small, and ps is not too large. The dependence on the 
concentration pa of fuel molecule for a fixed pb in this 
expression agrees with the one known as the Michaelis- 
Menten equation, and such a dependence of the rota- 
tion and reaction rates on the ATP concentration was 
observed for the Fi motor. 4 ' 27 The dependence on ps 
predicted in cq. (7) may be observed for the Fi motor as 
the dependence on the ADP concentration. 

Examples of the dependences of v and r on L arc 
shown in Fig. 2. Both v and r decrease monotonically 
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Fig. 3. The dependence of the efficiencies of motor (the upper left 
portion) and of generator (the lower right portion) on L and A/i 
is shown as a contour plot for pg = 1. The other parameters are 
the same as the ones in Fig. 2. In the hatched region the engine 
works as neither a motor nor a generator (it is useless). 



with increasing L. The rotation rate v becomes zero at 
a certain value Lq of L, and the reaction rate r becomes 
zero at a somewhat larger value L\ (see the upper inset of 
Fig. 2). Thus, the engine works as a motor for < L < Lq 
and as a generator for L > L\, whereas it wastes both the 
chemical and mechanical energies for Lq < L < L\. It is 
not difficult to see that Lq < Ap/2ir < L\ from property 
(iii) mentioned above, and that both L and L\ tends to 
lS.pl In (therefore the useless interval of L vanishes) in 
the tight coupling limit [exp(— ^o/k-gT) — > 0] according 
to property (iv). 

In the examples shown in Fig. 2, the mechanochemical 
coupling is almost tight {y ~ r) for small load. However, 
the difference between v and r becomes apparent (the 
coupling becomes looser) as L is increased. This is be- 
cause the effective barrier height $ — Wo — L(9b — 9c) 
of potential Vb seen from the location 9 = 9b decreases 
and hence the leftward dashed passes in Fig. 1 are taken 
more frequently for larger L. 

The efficiency n of energy transduction is defined as 
n = 2-KvL/r/S.p for motor and n = r/S.p/2irvL for gener- 
ator. 18 These expressions may be written as 



v = xm, 



(8) 



where \ — v l r { r / v ) IS the "tightness" of mechanochem- 
ical coupling, and 770 = 2ttL/Ap (Ap/2nL) is the ef- 
ficiency in the tight coupling limit for motor (genera- 
tor). Note that larger \ indicates tighter mechanochem- 
ical coupling, and we have \ = 1 in the tight coupling 
limit. 

An example of 77 as a function of L is depicted in the 
lower inset of Fig. 2, where we observe that n for motor 
(generator) has a maximum near the "stall" load Lq (Li). 
In the tight coupling limit, the maximum efficiency of 
rjo = 1 is achieved at L — Ap/2ir =p (the minus sign 
for motor and the plus sign for generator). In the case 
of loose coupling, the maximum value of 77 tends to be 
larger for larger \- Since the tightness decreases with 
increasing L as explained above, the maximum of n for 
motor is larger than that for generator in this example. 

The dependence of n on L and Ap for a particular 
choice of ps is shown as a contour plot in Fig. 3. In 
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of various results presented here, and its extensions to 
situations other than the rotary motor will be reported 
in future publications. 
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Fig. 4. The maximum efficiencies r? m (ps) of motor and of gener- 
ator for given ps are plotted against pB for different choices of 
<E>0. The other parameters are the same as the ones in Fig. 2. 



this example, the largest efficiencies of motor and gener- 
ator are achieved in the condition Afi ~ 2irL ~ 20fceT 
far from equilibrium (Afi = L = 0); a similar observa- 
tion was made for related models by Parmeggiani et a/. 18 
Note that, in the tight coupling limit, we have the max- 
imum value of r)o = 1 on the diagonal line Ap = lixL. 

We have obtained qualitatively similar patterns of con- 
tour lines to the one shown in Fig. 3 for other choices of 
ps, although the locations and the hights of the peaks 
change as ps is varied. Let rf m (pB) be the maximum 
value of 7] obtained by adjusting L and Ap for a given ps- 
Figure 4 shows Tjui cLS £L function of pb for three choices 
of the barrier height $o of potential Vq. It is noted that 
r) m increases with ps for both motor and generator, and 
saturates to a certain value — the upper limit of efficiency 
achieved by the engine (characterize by Vq(9), Vi(9), ka, 
kb and Do). 

The dependence of rj m (pB) on ps shown in Fig. 4 
may be understood qualitatively as follows. Remem- 
ber that larger rf m is expected for larger tightness \ °f 
mechanochemical coupling. As the the concentration ps 
of molecule B is increased, the binding of molecule B 
(transition from state to state 1 at 9 = 6b) occurs 
more frequently, which in effect leads to the decrease in 
the chance of taking the leftward dashed passes in Fig. 1 . 
Therefore, the tightness \ an d hence r/ m will increase 
with ps, which is consistent with what we see in Fig. 4. 
In the case of motor, the dissociation of molecule B oc- 
curs more frequently than the binding, and the former 
is not affected by ps, while the binding of molecule B is 
essential for generator. This explains the stronger depen- 
dence of rf m on ps for generator than for motor observed 
in Fig. 4. 

It may be worth remarking that the present model 
may be viewed as a motor driven by ion-flow across a 
membrane: 34 ' 35 the transition at 9 a can be viewed as 
ion exchange with the outside of the membrane, and the 
one at 0b as ion exchange with the inside; in this case 
only one chemical species is involved. 

In summary, we have proposed a minimal model for 
molecular chemical engines that properly takes into ac- 
count the effects of fuel and waste molecules. The model 
is simple enough to work out various properties of the 
engine such as the efficiency of energy transduction. The 
detailed analysis of the model, including the derivations 
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